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We study the chiral phase transition at finite temperature in the linear sigma model by em- 
ploying a self-consistent Hartree approximation. This approximation is introduced by imposing 
self-consistency conditions on the effective meson mass equations which are derived from the finite 
temperature one-loop effective potential. It is shown that in the limit of vanishing pion mass, namely 
when the chiral symmetry is exact, the phase transition becomes a weak first order accompanying 
a gap in the order parameter as a function of temperature. This is caused by the long range fluc- 
tuations of meson fields whose effective masses become small in the transition region, ft is shown, 
however, that with an explicit chiral symmetry breaking term in the Lagrangian which generates 
qs^ . the realistic finite pion mass the transition is smoothed out irrespective of the choice of coupling 

strength. 
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I. INTRODUCTION 



Chiral symmetry plays a vital role in low energy phenomenology of pion-pion and pion-nucleon interactions [0. 
The symmetry is manifest in the underlying QCD Lagrangian in the limit of vanishing quark mass, but it must be 
broken spontaneously in the QCD vacuum. Historically, the profound implication of this phenomenon to the very 
existence of the pion was first conjectured by Nambu long before the advent of QCD, in analogy to the BCS 
y—( • theory of superconductivity, and it led Goldstone (^] to explore more general consequences of the broken symmetry 
t-H \ discovering his celebrated theorem. 

\& . There are many reasons to believe that this hidden chiral symmetry becomes manifest at high enough temperatures. 

One can take an analogy to the spin systems in statistical mechanics, like the Ising or Heisenberg models, where the 
spontaneous magnetization disappears as the temperature of the system is raised above the critical temperature: 
this is a prototype of the phenomena generally known as the order-disorder phase transition j|. Field theoretical 
I analog of the order-disorder phase transition at finite temperature was first investigated by Kirzhnits and Linde 
in the context of electro-weak interaction in cosmological setting, and then a systematic method of computation was 
developed by Dolan and Jackiw Q and others |j. 

A great amount of works have appeared in the literature on the problem of the chiral symmetry restoration at high 
temperatures and at high baryon density. Lee and Wick [|| conjectured the possible existence of "abnormal state" 
of dense cold nuclear matter where the chiral symmetry is restored. Baym and Grinstein |j] studied the chiral phase 
transition at finite temperature by applying the approximation methods developed earlier in many-body theory. These 
early expositions of the problem used the linear a model [ [To| which is suited for studying the role of the quantum 
fluctuation since the model is renormalizabile at least at the level of perturbative computation . Non-perturbative 
aspects of QCD at finite temperature have been studied more recently by the Monte Carlo method for the Euclidean 
path integral formulations of QCD on finite size lattice; such studies have shown that the chiral transition is abrupt 
and closely related to the quark-liberating deconfining transition 113] . There are many other works based on various 
phenomenological chiral models such as the four-fermion (quark) interaction model of Nambu-Jona-Lasinio (NJL) 
type (^|l|; the non-linear sigma model of Weinberg 14pl| which may be considered as an effective theory of QCD 



at low energies. 

Recent renewed interests in the problem are partially due to the suggestion by Bjorken |16| and others |l7j] that 
there may be some observable consequences of the chiral symmetry breaking transition if it occurs in the course of a 
highly relativistic collision of hadrons or nuclei; a chiral condensate may occasionally grow in the wrong direction as 
a fluctuation during the rapid cooling of the dense hadronic matter (quark-gluon plasma) produced in such a collision 
and this may result in a significant fluctuation in the number ratio of charged and neutral pions. 
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In this paper we re-examine the chiral phase transition at finite temperatures in terms of the linear sigma model. 
Our purpose is to describe the phase transition in static equilibrium in a simple self-consistent approximation, which 
is well-known in the many-body theory as the Hartree approximation. We will not address more difficult problem of 
describing the dynamics of the phase transition jl8],[l9| in the present work. 

We adopt the linear sigma model since its simplicity to represent the chiral symmetry is most suited for our purpose 
of investigating the symmetry aspect of the phase transition. The model has an advantage of being renormalizable at 
least in the perturbative sense. The Lagrangian density is given by [E0| 



C = C s +C sb (1) 



with 



1 777 

C s = #7 • 8 + g{a + ztt • f 7 5 )]^ + -({da) 2 + (cW) 2 ) - — (a 2 + t? 2 ) 

-±(<r 2 +n 2 ) 2 , (2) 
C sb = -ea (3) 

where ip, a, and 7r represent the nucleon, sigma, and pion fields, respectively. The term C s is symmetric and 
invariant under an SU(2)l x SU(2)r chiral group. The right and left combination ip R = ^(l + r y 5 )ip, ip L = — 7 5 )V> 
transform respectively according to the representation (1/2, 0) and (0, 1/2) while the sets (a, 7?) belong to the (1/2,1/2) 
representation. C s b is the symmetry breaking term. Two Noether's currents associated with (|J), namely the vector 
current and the axial vector current, are given by 

-» T 

An = -i/j + od^ix - nd^a 

respectively. The equations of motion for the fields derived from the Lagrangian density (|J) give the PCAC relations 

= evf. (4) 

In the following investigation, we concentrate on the finite temperature behavior of the theory at zero net baryon 
density or zero chemical potential for baryonic charge and we will focus on the meson sector of the Lagrangian (||) 
which is 0(4) invariant under the rotation of the meson multiplet (a, tt). This may be justified at low temperature 
where the thermal creation of baryon-antibaryon pairs is suppressed due to the large baryon mass. This may not be 
the case, however, in the transition region where the effective baryon mass becomes small pi] . Inclusion of baryonic 
fluctuation is straightforward but we omit it here for the consistency of our approximation scheme (see Appendix A). 

All informations concerning the equilibrium properties of the system are contained in the thermal (imaginary- 
time) Green's functions defined by the ensemble average of the "time ordered" products of the quantum fields in the 
imaginary-time Heisenberg representation, cf>(x) = 4>(r,r) = e rH 4>(v)e~ TH |22] ]: 

n {m), \-/rrli \ li w Tr [^ PH T T ^{xi)---i{x m )] 

'(Xi,- ■ ■ ,x m ) = {T T 4>(x{) ■■■<t>(x m )) = (5) 

Ir e^P" 

where H is the Hamiltonian of the system and (3 = 1/fceT is the inverse temperature. The time ordering operator 
T T rearranges the field operator with the smallest value of r (0 < r < S) at the right of the sequence. The thermal 
Green's function can be expressed in the Euclidean path integral form p^pl[ |: 

G { ™\x u ---,x m ) = Zp 1 J V[<j ) ]<j ) (x 1 )---<j ) (x m ) e Mj^xC{^d<P)] (6) 

where £((f), d(f>) is the Lagrangian density of the system in the Euclidean metric, obtained by the substitution, 
d(f>(x)/dt — > id4>(x)/dr, and 

Z p = Tr e- 13 " = J V[<t>] exp[ J d A xC{(f), d<j))] (7) 
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is the partition function of the system. Here we have introduced a short hand notation for the Euclidean space-time 
integral: d 4 x = J d 3 r dr. In the above expression the functional integration over the classical field variable 
<p(r, t) defined in the range < r < (3 is performed with periodic boundary condition: <p(r, 0) = <p(r, ($). 

The real-time Green's functions, which describe the propagation of an external disturbance introduced in the 
system in equilibrium, may be obtained from the imaginary-time Green's function by the analytic continuation r — * it 
p2J. Therefore, the thermal Green's function also contains some dynamical properties of the system (for weak 
perturbations) such as the effective mass of the excitations which we wish to study in this paper. The well-known 
advantage of using the thermal Green's function is that it allows one to use systematic perturbation series expansion 
and evaluate each term by the method of the Feynman graphs. 

In the next section we review the effective potential and its loop expansion at finite temperature. We will apply 
the method to calculate the effective mass of the mesonic excitations at finite temperature in section 3 and derive the 
Hartree approximation. The resultant gap equations, or the self-consistency conditions on the mass parameters for 
sigma mesons and pion, contains ultraviolet divergent integrals associated with the vacuum fluctuation and therefore 
requires renormalization. We encounter a difficulty in imposing renormalization conditions to absorb the infinities in 
the redefinitions of the mass parameter and the coupling constant as in the perturbative renormalization scheme: the 
only consistent renormalization conditions in our nonperturbative approximation scheme give symmetric solutions. 
This difficulty led us to adopt the phenomonclogical approach to ignore the vacuum polarization effects and to retain 
only thermal fluctuation in calculating the loop integral in the gap equation. 

In section 4 we examine the numerical solutions of the resultant gap equation. We will show that they exhibit 
a nontrivial hysteresis behavior characteristic of the first order phase transition in the exact chiral limit. It will be 
shown that this hysteresis disappears when we introduce explicit symmetry breaking term consistent with PCAC 
irrespective of the choice of coupling. 

A short summary and conclusions are given in section 5. Two appendices contain the discussion on the inclusion 
of the baryon fluctuation and the high temperature expansion of the loop integral. Throughout the paper we use the 
natural unit: c = h = 1. 

II. THE EFFECTIVE POTENTIAL AT FINITE TEMPERATURE 

In this section we review the effective potential formalism and the loop expansion at finite temperature. The 
finite-temperature effective potential for the linear a model is given up to the one- loop order. 

A. Some general formalism 

The thermal Green's functions (4) can be obtained from the generating functional W^[J]: 

e w P = Zp{J] ee (T T exp[ / d 4 x^(x) J(x)]) 

= Zp 1 J V[<t>] exp[ J d 4 x (C(<f>, d4>) + <j>J)] (8) 

by the functional differentiation with respect to the "external source" J(x). For example, the thermal average of the 
field operator <f>(x) is given by 



<*.» - 5m ' 1 



5J{x) 



(9) 

j=o 



and the other thermal Green's functions can be extracted from the following Taylor series by proper functional 
differentiations: 

Wp{J) = Y,~, [ d 4 x 1 ---d 4 x m G ( J H \x 1 ,---,x rn )J(x 1 )--- J(x m ). (10) 
„=i n - h 

The effective action r^[0] is defined by the (functional) Legendre transform of the generating functional: 

VS\ = Wp[J] - [ d 4 x4>{x)J{x) (11) 

h 
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where J(x) on the right hand side is considered to be a functional of fax) determined by fax) = SWp[J]/SJ(x). It 
follows then that 

and therefore the thermal expectation value of the field (8) corresponds to a stationary point of the effective action, 
while the "curvatures" of the effective action at such local minimum are related to the two-point thermal Green's 
function: 



S(f>(x)5(f>(x) 



SJ(x) 



J=0 *m 



= -Gf{x,y)-\ (13) 

j=o 



In the following analysis, we assume that the system is homogeneous and therefore translationally invariant as in the 
vacuum. In this case the expectation value of the field becomes independent of x, and the two-point thermal Green's 

(2) 

function Gp (x,y) becomes a function only of the difference of two space-time arguments; its Fourier transform is 
therefore given by 

=.(2) 



G^\k)= / d*xe^G^{x,y) (14) 

where k = (ui n , k) and uj n — 2irnj3~ l with integer n is the discrete Matsubara frequency. 
The effective potential Vp is defined by the effective action with the constant fax): 

Vp{<t> ) = -Tp[4>]/n for fax) = fa (15) 

where f2 = J_ d A x = (3 J d 3 r is the four-volume of the Euclidean space-time. The thermal average of the field is now 
determined by the stationary condition of the effective potential: 

5M« =0 (16) 
ofa 

Since the effective potential corresponds to a part of Tp[(j>] which contains only k = mode of fa = Jp d 4 x(j)(x)e lkx , 
its second derivative with respect to fa gives the Fourier component of two point Green's function evaluated at zero 
external momenta: 

^S0W-a?(p)-'. a-) 

The four point vertex function Xp(ki, k2,k%, k±) at zero external momenta k\ — ki — k^ — k^ = is also determined 
by the fourth derivative of the effective potential: 

^#^ = A,(0, 0,0,0). (18) 

The path integral expression (||) of the generating functional W^[J], together with the definition of the effective 
action (pi]), allows one to derive a loop expansion of the effective potential Vp p^,p[: 

Vpifa) = Vo(fa) + Vf(fa) + AV?(fa), (19) 

where the first term Vo(fa), the tree approximation, is the potential term in the classical Lagrangian. This term is 
independent of the temperature. Temperature dependent one-loop term Vf and higher loop corrections AV 13 are 
obtained by the following procedure: We decompose the original Lagrangian by shifting the field, fax) = fa + 4>'(x): 

£[0] = C [4>'\ fa] + £'iW\ fa] — Vq + linear terms + constant, (20) 

where C' [4>',fa] is a bilinear form of the shifted fields 4>'{x) and the new interaction term C'j contains higher order 
products of <p'(x). The one-loop contribution V\ to the effective potential is then given by 
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e -nv M) = / v[4>l] exp[ / d 4 xW . (21) 



Since this functional integral Gaussian it is easily performed: 



Wo) = -5 J^d 4 k\nD(k,M, (22) 

where we have introduced a short hand notation for the integrals over the momenta and the sum over the discrete 
Matsubara frequencies: 

and the "free" propagator D(k, fa) is defined by 

/ d A xC'^fa,fa] = I ^k\D-\k,4> )$i. (23) 

J/3 J/3 1 

The multi-loop corrections AV^ is given by 

-oav(^) = /P[0 / ]exp[^^(£-(^;0 o )+r / (^;0 o ))] 

which can be interpreted diagrammatically as the sum of all single particle irreducible connected graphs involving the 
vertex functions given by the shifted interaction Lagrangian C'j with no external leg and internal lines are associated 
with the "free" propagator D (k : (j) ). 

In passing we note that the effective potential can be physically interpreted as the Hclmholtz free energy density 
of the system in the presence of "self- magnetization" fa: 

Vp(fa)=a(T,fa)=g(T,J)-faJ, (25) 

5 (T ) J) = -^lnTre^^-/^ (r)J ], (26) 

where g(T, J) corresponds to the Gibbs free energy density of the system with the "external field" J. From the 
thermodynamic relations, it is related to the pressure P(T, fa) of the system by 

Vp(fa) = -P(T,fa). (27) 



B. Effective potential for the linear a model 



The leading terms of the loop expansion of the finite temperature effective potential for the O(N) a model was 
computed by Dolan and Jackiw || . If we ignore the baryon sector in the a model Lagrangian, the meson sector of the 
Lagrangian has the 0(4) symmetry with respect to the rotation of the meson fields (f> a (x) = (<j(x), tti(x), ^2(2;), ^3(2;)) 
in the limit of vanishing pion mass. This symmetry is directly reflected in the effective potential Vp(4> a ) at all 
temperatures and the spontaneous symmetry breaking is signaled by the appearance of the minima of Vp((f> a ) at 
non- vanishing cf> a . 

The leading term of the loop expansion of the effective potential is given for the linear a model by 

V (4> a ) = ±m 2 4> 2 + ±\^-ea. (28) 

This term is independent of temperature and contains the symmetry breaking term proportional to a. As we shall 
see below this is the only term which depends on the symmetry-breaking term in the Lagrangian C s b : all other terms 
due to quantum fluctuations therefore should depend only on cp 2 — 4> a 4>a- 

To compute the one- loop contribution V{ ( and higher loop corrections AF' 3 ), we shift each component of the fields 
by a constant amount <f> a : Q (x) = 4>' a (x) + <pa and decompose the Lagrangian density: 
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C[<j) a ] = jCo(0q(x), <f) a ) + Cj(<p' a (x), cj) a ) + linear terms + constant, (29) 
where the shifted "free Lagrangian" Co(4>' a (x), <f> a ), which is bilinear in the shifted fields cf>' a , is given by 

Co(<f> a (x), 4> a ) = x {d4> a d<t) a - <p a m 2 ab <j>b) (30) 



with the mass square matrix given by 



m lb = I™ 2 + 7rA0 2 ](5 afa + iA0 a & , (31) 
o 3 



where <j) 2 = ^ a 4> 2 an d hereafter all repeated field subscript imply the sum if not otherwise indicated. The shifted 
interaction Lagrangian becomes 

Ci(</>a(x), 4> a ) = ~\\^ a 4>a{x)4>' 2 {x) ~ ^M\x). (32) 

6 4! 

The shifted free Lagrangian Co determines the free propagator at finite temperature: 

d i xL o (<t> a (x),0 a ) = l [ d*kMk)\p-HWa)UMk) ( 33 ) 

1 J (3 

where [D~ 1 (k, 4> a )]ab = k 2 S a b + m 2 ab with k 2 = w„ + k 2 . The diagonalization of this matrix can be performed by 
rotating the fields to cj)' a = O a b4>b = (4>, 0, 0, 0): 

fc 2 + m 2 



k 2 + m 2 



where 
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m\ = m 2 + ]^\<j> 2 (35) 

TO 2 . = TO 2 + — \<j) 



and wc find 



G 



D ab {k;^ a ) = + _ 2 S ab - -=z- ■ (36) 



k 2 + TO 2 2 ru -r iit 2 

The one-loop contribution to the effective potential is given by 



tf(<« = J / ^fclndctD^;^)] 

1 Jf3 

= -^ _1 E / (0Kfc 2 +-i) + 31n(fc 2 + m 2 )] (37) 

where we have indicated the sum over the Matsubara frequency ko — uj n — 27m/3 -1 explicitly. This sum may be 
performed by the following formulae: 

E (2n7r) 2 + x 2 = 2 + 1 ^ 

n 

^ In [(2mr) 2 + x 2 ] = | + In 1 1 - e~ x \ + cq (39) 

n ~ 

where the formula (|3^) can be obtained by the method of contour integration |^,D while ( p9| ) is derived from ([38j) by 
integration over x. (The constant of the integral Co = \ J2 n ^o l n (2?w) 2 is actually infinite. ) The result is summarized 

as 



G 



(40) 



where 



d 3 k 
(2^)3 



[e k (mi) + 3e k (m 2 )] 



d 3 k 
J2nf 



ln(l 



-f}ek(mi)\ 



31n((l 



-,8eic(m2) 



(41) 
(42) 



with £ic(m) = + to 2 . The first term V^ is due to the vacuum fluctuation; this term is divergent but is temperature 
independent so that it may be removed by some renormalization procedure at zero temperature. The temperature- 
dependent term Vf can be rewritten by integration by parts: 



d 3 k 

(2^)3 



3e k (TOi) 



/B(£kK)) 



k 2 



3£k(TO 2 ) 



/B(£k(m 2 )) 



where 



ME) 



1 



(43) 



(44) 



is the bose distribution function. Recalling that the effective potential is equal to the pressure of the system with 
inverted sign, this term may be interpreted physically as due to the pressure exerted by the ideal gas of a boson with 
mass toi and another kind of boson m 2 with three-fold degeneracy. 

The multi-loop corrections AV@ may be computed from the formula ([m]) by using the shifted Lagrangian densities 
C'q and C'j. This term may be interpreted as the contribution of the interaction among these quasi-particles (defined 
by C ) due to the residual shifted interaction given by C\. We note that both C' and C\ do not depend on the 
symmetry breaking piece C s b of the original Lagrangian density. It therefore follows that AV@ should only depend 
on dr. 



III. SELF-CONSISTENT HARTREE APPROXIMATION 



In this section we derive a self-consistent approximation scheme for the computation of the thermal Green's function. 
To motivate the method we first examine the effective masses of the mesonic excitations at finite temperature using 
the finite temperature effective potential in one-loop approximation. The result is used to introduce a self-consistent 
Hartree approximation at finite temperature. We show a difficulty in the renormalization procedure and this leads us 
to take a phenomenological approach to redefine our Hartree self-consistency conditions. 



A. Temperature-dependent effective mass in one-loop approximation 

The effective potential V(4> a ) depends only on (jr in the absence of symmetry breaking term C s b in the Lagrangian. 
When wc include the symmetry breaking term C s b in the Lagrangian, the only term which breaks this symmetry is the 
last term in Vb . It then follows that the stationary conditions with respect to the variations of <p a become degenerate 
except for the a field: 

/ i »t^/3 \ 

e = 0, (45) 

= 0. (46) 

It is clear that the potential minimum always appear at ffj = 0. In the following, we therefore assume ffj = and 
consider only the first condition ( |45| ) which determines the potential minimum at a — <7q = y^o- 

Having determined the equilibrium conditions for the condensate amplitude, we now compute the second derivatives 
of the effective potential at the potential minimum; this determines the inverse two-point thermal Green's functions 
in one-loop approximation at k — 0, which correspond to the mass square of the excitation, Gp a ,^(k = 0) = Af 2 ^: 
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d 2 vf 

a— ctq ,7f— 



Ml = m? + I ) + - (47) 



M OT 2 = mi+[^n +-.. (48) 



<9tt 2 



a— ao ,7r— 



where the first terms came from the differentiations of the tree effective potential Vq and the explicit form of the 
contributions from the one-loop term ( |37j ) are given by 



da 2 I J \ k 2 + m\ k 2 + : 



v 2 -2 ' 1 1 



- A ^y/Ho?w + 3(i?+^j' (49) 

) = A / d 4 k ( -„ 1 . + ,„ 1 . (50) 

\ / CT — (Tq ,7T — 



Since 



daditi 



| =0 (51) 

(T— (7Q ,7f— 



there is no a and 7r mixing even at finite temperature. 
We note that since 



' hT| ' _ ^A / I — ; : -r + ttt - "j ) (W. 



the stationary condition (^) for the cr field gives the following relation for the pion mass 

M 2 ^{(i) = e (53) 

in the one-loop approximation. In the vacuum (/? = oo) the condensate amplitude cto is equal to the pion decay 
constant f n defined by the transition matrix element of the axial current (0\Aj l (x)\n k (q)) = i5jkq f _ l f 7r e~' lq ' x since 
A 1 ^ = ad^TXi — Ttid^a in our model. Therefore we have 

m 2 n U = e (54) 

in the vacuum. The above relation guarantees that if we have exact 0(4) symmetry in the Lagrangian (e = 0), the 
pion becomes massless when a ^ 0, consistent with the Goldstone theorem. 



B. Self-consistency conditions for the Hartree approximation 



The equations ([17]) and ( pq ) are nothing but the Schwinger-Dyson equations for the two-point thermal Green's 
functions evaluated at zero momentum: 

G^ A (fc = 0) = D~^(k = 0) + n M „(k = 0) (55) 

where IIg )0 . /„.(&) is the self-energy of a and pions with Euclidean momentum k = (u;„,k) respectively. In the one- 
loop approximation these self-energy terms are evaluated for the one-loop Feynman graphs shown in Fig. ([l]) using 
Dp t<T / n (k) for the internal propagators. 
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(C) (d) 

FIG. 1. One-loop self-energy diagrams for a, (a) and (b), and for pions, (c) and (d). Solid line is the free sigma meson 
propagator and the dashed line is the free pion propagator. 

This observation suggests j(| that a partial sum of the infinite series of the loop expansion may be made by imposing 
the self-consistency condition on the one- loop result replacing the propagators D a /^(k) for the intermal lines of the 
one-loop Feynmann diagrams by the Hartree propagator of the form 

Gl^{k) = \e + Ml / X X . (56) 

We define our "Hartree approximation" by imposing self-consistency conditions for the mass parameters M a /^ intro- 
duced in this propagator: masses 

Ml = m 2 + ^ 2 + n^[M a ,M^, (57) 

Ml = m 2 + ^\4? + n^(M CT ,M w ) (58) 

where two self-energy terms are given by 

nl a (M a ,M„) = U^M,, M w ) = A [^(M 2 ) + ^(M 2 .)] (59) 

with 

^ M2 ^j p ,/ k WTW (60) 

where the subscript B is just to remind that the sum is taken over the boson Matsubara frequency These are the 
generalization of the Dolan-Jackiw "gap equations" for multi-component fields in the presence of non-vanishing con- 
densate ap = 0. This approximation is called "modified Hartree approximation" by Baym and Grinstein. Successive 
substitutions of the right hand sides of the equations into the arguments of the self-energy terms generate an infinite 
series of "superdaisy diagrams" . 
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(a) 



71 



(b) 

FIG. 2. Schwinger-Dyson equations in the self-consistent Hartree approximation for thermal sigma meson propagator (a) 
and for thermal pion propagator (b). 

We note that the direct substitution of the internal propagators in the one-loop self-energy for the a field would 
generate an extra term for the a self-energy Iljj a , 



where 



: Ml) + : Ml) 



1 



^(fc : M 2 ) _ d 4 k' \ k a +M 2 ){{k , + k) 2 +M 2y 



(61) 



(62) 



We have excluded this term since it has non-trivial k dependence (or "dispersion" ) and therefore invalidates our simple 
ansatz ( |56|) for the Hartree propagator. 




FIG. 3. The diagrams omitted in our approximation. 

Our prescription for the self-consistent Hartree approximation is not complete yet since we have not specified the 
condition for determining the condensate amplitude a which appears in the self-consistency condition for ^l 1 p~^°° p - 
Our guiding principle here is the Goldstone theorem. We replace the one-loop stationary condition for the a field by 



m 2 + -A ( 7 2 +nH 7r (Af ff ,M /3 ) 



a = e. 



This equation, together with (EH), implies 



Mia = e 



(63) 



(64) 



so that it guarantees that the pion becomes massless Goldstone mode in the limit of e — > 0. 
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C. Renormalization 



The self-energy terms in one-loop approximation and in self-consistent Hartree approximation contain (quadratic) 
divergences which all originate from the divergent integral of <&@(M). These divergence must be removed carefully 
by the renormalization procedure due to the nonlinear structure of the equations. It turns out, however, that such 
non-perturbative renormalization is possible only when there is no symmetry breaking a = as we shall see below. 

We first isolate the divergence in the self-energy terms by performing the Matsubara frequency sum using the 
formula Q: We find 



(2tt)3 2e fc 



•$/j(M) 



(65) 



where 



d 3 k 1 



/B(£k). 



(66) 



All divergences in the self-energy terms arise from the phase space integral of the first temperature independent term. 
The phase space integral of $>p(M) converges due to the ultraviolet cut-off by the distribution function /e(efc) and 
this term vanishes at zero temperature. To regularize the divergent integral we introduce the ultraviolet momentum 
cut-off A in the phase space integral: 



rf J fc 1 



7 1 (A)-M 2 / 2 (A,/i) 



M 2 , M 2 
■ m — 7t 



16tt 2 



0(A~ 



(67) 



where 



/ 2 (A )M ) 



A2_ 

87T 2 ' 
1 

16tt 2 



[ln(4A 2 /M 2 ) - 1] 



(68) 
(69) 



and /i is an arbitrary constant which plays the role of a renormalization scale. 

We seek the renormalization conditions for the bare parameters (m, A) so that the resultant equations are written 
only in terms of the renormalized parameters {tur, Xr) and contain only finite terms. For this purpose it is convenient 
to rewrite the equations (^) and (|58| ) in a more symmetric form. We first add the two equations: 



Mi 



Ml = 2m 2 



= 2m 2 



^Xa 2 



2A [$f}(M a ) + ^(M,)] 



+2A 



16tt 2 



■hi 



2 

3' 

Ml 



4A/ X (A) + -Xa 2 - 2AJ 2 (A, M )(M 2 + M 2 

' M 2 



M 2 M 2 
" hi 

V 2 



16tt 2 



2A 



$^(M CT ) + $ /3 (M 7r ) 



(70) 



where we have used (B7]) in deriving the second line. This equation may be rewritten as 



M 2 + M 2 = 



-2m 2 R 



2 \ a 2 



2Xj 



Ml M 2 a 
in — — 



16tt 2 



M 2 M 2 
- ln^f+^(M CT )+^(M 7r ) 



16tt 2 



(71) 



by imposing the following renormalization conditions pi 





2 

m 


Xr 


= X 


1 


1 


Xr~ 


= A 4 



(72) 
(73) 



On the other hand, however, the subtraction of (58) from (|57|) gives another relation 



Ml - Ml = -Aa 2 



(74) 
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which contains the bare coupling A. In the limit of A 
we find 



oo, the bare coupling goes A 



and in this case the renormalized self-consistency condition (71) becomes simply 



M 2 = -ml 



2\r 



M 2 M 2 ■ 

In — y + ®p (M) 



16?r 



0_ for any finite A#, hence 
(75) 

(76) 



for M = M a = M n . Unfortunately, this renormalization scheme is possible only when the symmetry is not broken 
ct = and e = and does not apply for the broken-symmetry phase. 

Similar difficulty was noted in for the finite temperature self-consistent approximation similar to ours. 

Having seen the difficulty in finding the decent renormalization conditions for our self-consistency conditions, we take 
an phenomenological approach to the problem: we replace the divergent integral ^p{M 2 ) by its finite temperature- 
dependent piece ^(M 2 ) simply discarding the divergent integral associated with the "vacuum loops". We therefore 
impose 



Mi 



Mi 



\M 2 + A 
+ A 

6 



5/3 (AO 



m 2 + -\<f + A \ $p(M£) + 9f,(Mi) 

for the condition to determine the amplitude of the condensate. 

The physical rationale for this procedure is the observation that each term in the self-energy term, e. g. 

d 3 k 



for the temperature-dependent effective masses and 

1 



AS^Af ) = A 



1 



-/ B (£k(M)) 



(77) 
(78) 

(79) 
(80) 



(2tt) 3 2e k (M)- 

is just due to the forward scattering of the meson by other mesons of effective mass M present in the system. Since 
the a model should be considered as an effective theory of the underlying microscopy theory of the strong interaction, 
namely QCD, it would be physically sensible to only take into account the effect of physical excitations in the system 
as these equations actually do. This approximation also corresponds to neglecting the "many-body" interactions 
among the thermally excited mesons due to the vacuum loops. 



IV. NUMERICAL ANALYSIS OF THE GAP EQUATIONS 

In this section we present numerical solutions of the gap equations derived in the previous section. It will be shown 
first that in the exact chiral limit the solutions exhibit hysteresis behavior characteristic of first order phase transition. 
We will show, however, that the introduction of the finite symmetry breaking term washes away this hysteresis for 
any choices of model parameters in accordance with the physical conditions. 

In constructing numerical solutions of the gap equations (|77|), ( f78| ) and (|79|), we first note that from ([77j) and ( |78| ) 
the effective sigma mass is related to the effective pion mass by 

Ml = Ml + ±\a 2 (81) 

while ( ]78| ) and ((79|) imply that the effective pion mass is related to the strength of the explicit symmetry breaking 
term in the Lagrangian density by 

Mia = e. (82) 
What remains to be done is to solve one of the self-consistency condition, say for the effective sigma mass, 

Ml = m 2 + X -M 2 + A [^(M 2 ) + ^(M 2 )] , (83) 

by inserting @ and @. 

Although the physical value of e is actually given in terms of the physical pion mass and the pion decay constant 
by the zero temperature relation (e = to 2 ./^), we consider e as an external variable and examine how the solutions of 
the gap equations depend on the strength of e. 
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A. Exact chiral limit (e = 0): 



We first examine the special case when e = 0. In this limiting case the chiral symmetry is the exact symmetry 
of the Lagrangian and according to the Goldstone theorem the pion becomes massless in the low temperature phase 
where the symmetry is spontaniously broken. Our approximation scheme in fact guarantees this condition by (|8T 
At zero temperature, the above relations (plj) and (p3) imply 



M 2 (T = 0) = mg = -2m 2 , 
a 2 {T = 0) = (g = -6m 2 /A. 



(84) 
(85) 



In the symmetry broken phase where a ^ 0, M n — and M 2 = ^Aer 2 , the self-consistency condition for the effective 
sigma meson mass becomes 



Ml = m 2 + ?M 2 



A 
2^2 



I m (M a ) + /i 2) (0) 



(86) 



where we have introduced dimensionless variable M a = M a (3, rh — m/3 and the dimensionless function (/i) is 
defined by 



/£ 2 V) 



x dx 



o \Jx 2 + y? e V x2 +^ 2 ± 1 



(87) 



The high temperature expansion of I±\p,) appropriate for small value of /i is given in the Appendix B. On the other 
hand, in the high temperature symmetric phase, where a = and M a — M n , the gap equations become degenerate 
to 



Ml =rh 2 + 2\ 7 ^I m (M a ) 



Using I±\o) = 7r 2 /6 the above two conditions are further simplified to 



6 



for the low temperature phase, and 



7T Z A 



12 



■k z A 



1 ~„ 

-Ml 



-m 2 + Mi 



(89) 



(90) 



for the high temperature symmetric phase. 

It is instructive to examine the solutions of the above equations graphically. In Fig. |] we plot the function 

f{x 2 ) = \l<?\x), 



(91) 



which appear on the left hand side of the above two equations, by the solid line; the right hand sides of (|89j) and 
are plotted by dashed curves. The intersection of two curves determines the solutions of the scaled gap equation, 
or (ficj), for the scaled sigma mass square M 2 . We note that f(z) is a concave, monotonically decreasing function of 
z = M 2 normalized as /(0) = 1, while the right hand sides of (p9|) and (90) are both linear in z. Since the value of 



-m 



is fixed by the amplitude of vacuum condensate at 



-m 



l/6A0g(> 0), the y-intercepts of the dashed curves 



(— 12?fi 2 /A — 1 for the low temperature phase; — 6m 2 /A for the high temperature phase) both move upward as the 
temperature is lowered, while the solid curve remains unchanged. 
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FIG. 4. Graphical construction of the solutions of the gap equations (gg) and (gfj): the solid curve depicts the left hand side 
of the equations and two dashed lines marked by LT and HT are for the right hand side of (^) and (po|), respectively, both as 
a function of M%. Crossing points of the solid curve and the dashed line give the solution of the corresponding gap equation. 

In the high temperature phase, two curves cross when — 6m 2 /A > 1. The condition — 6m 2 /A — 1 determines the 
"critical" temperature 

T = O (92) 

below which the real solution of ( |90| ) docs not exist. Note that To does not depend on the strength of the interaction A. 
When one approaches to Tq from high temperature side, the effective sigma mass decreases continuously and vanishes 
at T = Tq. On the other hand, at low temperatures below To, where — 6m 2 /A > 1, the dashed curve for the low 
temperature phase intersects with the solid one at one point giving a unique solution to the original gap equations. 
As the temperature increases, the effective sigma mass again decreases. However, it does not vanish at To in this case. 
Instead, as the temperature increases slightly above T , there appears another crossing point at small M 2 due to the 

concave shape of the function as plotted against ^ 2 . As the temperature is further increased, the two crossing 

points approach toward each other and they eventually annihilate at certain temperature Ti(> To) when the dashed 
line becomes tangent to the solid curve. 

We show in Fig. |5| the temperature dependence of the effective sigma mass and the amplitude of the sigma 
condensate, normalized by their zero temperature values given by ( j34| ) and (|8^) respectively. This plot is made by 
solving ( |89| ) and (|9(]) respectively for —rn 2 and making parametric plot of the temperature (oc l/V— m 2 ) vs the scaled 
sigma mass (oc M a j V—m 2 ). The scale of the temperature in this plot is set by the value of To which is equivalent 
to the amplitude of the vacuum condensate by (|9^) . The system exhibits a typical hysteresis behavior of first order 
phase transition in the exact chiral limit. We note that this hysteresis appears always irrespective to the value of A 
since the slope of the function f[z) diverges at z = 0. The magnitude of Ti depends on A explicitly, however. For 
larger A, the slope of the dashed lines in Fig. [| (b) becomes less steeper, thus T\ is larger. Since T is independent of 
A, this implies that the region of the hysteresis is larger for stronger coupling. We next examine how this behavior is 
modified when e is finite. 
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1.2 1.4 



FIG. 5. Temperature dependence of the solutions of the gap equations in the chiral limit (e = 0) for the scaled sigma mass 
(fj, a = Ma/mo) and the scaled sigma condensate (\ — ^/0o) as a function of t — T/Tq. In the high temperature phase, the 
effective pion mass becomes degenerate with the effective sigma mass while it vanishes in the low temperature symmetry broken 
phase. The scale of the temperature is set by To — </>o- 



B. For e 0: 



As we have noted earlier, e plays the role of externally applied magnetic field in the case of magnetic phase transition, 
therefore we expect that the transition will get smoother for non-vanishing e . 

For our numerical analysis with non- vanishing e, it is more convenient to use the following dimcnsionlcss variables: 



Ho = M a /m , 
fi„ = M n /m , 

X = &/<t>o, 
t = T/T = T/fo, 
e = e/(m%<f>o), 



(93) 
(94) 
(95) 
(96) 
(97) 



where m and 4> are the sigma mass and the sigma condensate at zero temperature in the limit e = as determined 



by (UJ) and (p5|), respectively. The equations (|8l]), (|82]) and (|||) then reduce to 



2 2,2 
= + X , 

Ax = e, 



/ 2) (v/A/3^A) +/^(v / A/3^A) 



(2), 



(98) 
(99) 

(100) 



respectively. Note that in this form the A dependence of the gap equation is absorbed into the A dependence of 

(2) 

the argument of the dimensionless function /_ . The solutions of these coupled nonlinear equations depends on two 
dimensionless parameters, A and e. In the limit e — > they coincide with the previous results. 

We show in Fig. |^ the e-dependence of the scaled order parameter \ plotted as a function of the scaled temperature 
t. It is seen that for small values of e the solutions exhibit hysteresis behavior (the back-bending shape of x)- As e 
increases, however, the curve stretches out gradually and for a large value of e, y b ecomes a monotonically decreasing 
function of temperature. From an inspection of the equations d98|), ( |99| ) and ( |100| ), we can see that the effect of e is 



expected to become significant when \ 



e/x in (| 



r(2) 



or yj \/3(in/t ~ 1 for the argument of the second I_ ' m 
(100); the former condition gives e ~ \ 3 anc l the latter e ~ 3t 2 x/A for the conditions that a significant modification is 
caused by e. This explains why for a relatively small value of e, say e = 0.001, we see in the plot that a considerable 
modification appears at small X, X < 0-1) an d that for larger A the effect is more significant. We found that the 
critical value of e for the disappearance of the trace of hysteresis is e c = 0.002 for A = 100. 
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FIG. 6. Temperature dependence of the scaled order parameter \ = °"/°"o- The curves are labeled by the value of e. Two 
figures correspond to two different choices of the coupling strength A. 



The physical value of e may be determined by adjusting the values of our model parameters (A, 



and e) so that 



they reproduce the physical pion mass (m n — 140MeV), and the pion decay constant / w = 93MeV and the sigma mass 
m a in the vacuum. In the linear sigma model f„ is identified as the amplitude of the vacuum condensate: <jq = f„. 
Then (|Sl|) and (p3) give constraints: 



m 



which implies 



and 



A = 3(m 2 - mD/fl 

— m 2 = (m a — 3m 7r )/2. 
2 > these two relations give a constraint A > 6m 2 // 2 



(101) 



(102) 



(103) 



We note that since — m > these two relations give a constraint A > 6m 2 // 2 = 13.6 on the allowed values for 
A. When we vary A from 50 to 500, m a varies from 400 MeV to 1200 MeV constrained by these relations. The 
corresponding values of e = e/mp^o are determined by using m 2 , = —2m 2 and (f>o = fn/xo where \o is the solution of 
the scaled gap equation at zero temperature. Some representative results are tabulated in Table [|. For large coupling 
A, the sigma mass increases as m^ ~ mo ~ yj A/3/,r and the relative strength of the symmetry breaking term decreases 
as e = (^ T /™o) 2 (/ir/ ( / , o) ~ 3A _1 (m 7r // w ) 2 . We show in Fig. ^ the results for the temperature dependence of the 
mass parameter for two different choices of the coupling strength A. 
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FIG. 7. Temperature dependence of the effective sigma mass (M CT ) and the effective pion mass (M w ) for e (=mlU). Two 
figures correspond to two different choices of the coupling strength A. 
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An interesting question is whether the hysterisis can survive for a sufficiently large value of A since e decreases for 
large A. This is not the case, however, since the critical value of the e also decreases for small A. One can see this point 

(2) C k / 

noting that the argument of the second I_ in (|100|) becomes independent of A at large A: ^X/i^/t ~ (m n / f n ) / \/xt 
as A — * oo. For any value of A the introduction of the symmetry breaking term e = m^f^ with physical value for m v 
and / w washes away the hysteresis behavior in the solution of the gap equatioon. 



V. SUMMARY 



In this paper we have studied the chiral phase transition at finite temperature using the meson sector of the linear 
sigma model. We formulated a self-consistent approximation scheme starting from the one-loop effective potential by 
imposing self-consistency conditions on the effective meson mass equation, keeping the meson self-energy independent 
of meson momentum. The resultant equations are just the Dolan-Jackiw gap equations for multi-component fields in 
the presence of non-vanishing meson condensate. In this approximation, the meson mass parameter in the thermal 
(imaginary time) meson propagator may be identified by analytic continuation as the mass of the real mesonic 
excitations introduced in the system in equilibrium although in more general cases such simple identification is not 
possible. This non-perturbative approximation scheme has, however, a difficulty in choosing proper renormalization 
conditions to eliminate the divergent loop integrals. We therefore adopted a phenomenological approach which just 
ignores the divergent vacuum loops. This procedure is technically equivalent to introduce a zero momentum cut-off 
in the vacuum loop integral. 

We showed that the solutions of the resultant gap equations does not reproduce usual second order phase transition. 
We found, instead, that they exhibit hysteresis behavior characteristic of the first order transition in the chiral limit. 
It is shown that this "premature transition" is caused by the long range fluctuation of the mesons fields whose effective 
masses become small in the transition region. We are however not able to calculate the transition temperature from 
the gap equations alone; to do so we need the information of the effective potential associated with the approximation 
we employed in this work. Such calculation may be performed by the method of the composite operator effective 
potential developped by Cornwall, Jackiw, and Tomboulis p7j . We have made a preliminary investigation in this 
direction |2^| . The result will be reported elsewhere. 

Inclusion of the symmetry breaking term generally tends to smooth sharp phase transition, but in general, the first 
order transition may survive if the symmetry breaking scale is small enough. Wc found, however, that this is not 
the case at least in our calculation: the physical pion mass, together with the PCAC relation for the axial current, 
gives a very strong constraint on the choice of the parameters of our model and the hysteresis behavior is smoothed 
out irrespective of the coupling strength. We found that smooth chiral transition takes place at T = 150 ~ 200MeV 
in our model, consistent with recent results from the state-of-art lattice QCD calculations at finite temperature jt2) . 
An advantage of the present approach based on the effective degrees of freedom, expressing crucial aspects of the 
symmetry behavior of the system, is that it is more tractable for investigating more difficult problem of the dynamics 
of the chiral transition in high energy nuclear collisions where we expect non-equilibrium aspects may play essential 
role J]~7| [l9|] . We also plan to investigate this problem in the future. 
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APPENDIX A: INCLUSION OF BARYON 



The original a model Lagrangian contains the baryon (nucleon) fields coupled to meson fields by Yukawa coupling. 
At low temperatures (T « m^), we may expect the contribution of the thermal baryon- antibaryon excitations is 
negligible. This may not be so however near the transition temperature since the baryon mass becomes smaller due 
to the reduction of the a condensate which dynamically generates the baryon mass. For example, in the lowest order 
of the Yukawa coupling g, 



17 



mjsr = gcr, 



(Al) 



(see Fig. baryon-self). We show in this appendix that the inclusion of the baryon Auction in this approximation leads 
to the results very similar to the Hartree approximation without baryons. However, we show also that it is not a 
consistent approximation for the meson propagators. 



X 



FIG. 8. The lowest order baryon self-energy diagram 

Inclusion of the baryonic excitations in the computation of the effective potential may be performed in the line 
similar to mesonic excitations (Sect. 2) with one modification that the functional integral over the Grassmann 
variables, which represent "classical fermion field", leads to the change of the sign and the removal of the factor 1/2 
in the one-loop contribution: 

^fbaryoB^) = / d*khl dsb G N (k) (A2) 

where 

G N (k) = 1 (A3) 
jk + M N 

is the thermal Green's function for baryons with = iu> n jo + %ki and the sum is taken over the fermion Matsubara 
frequency tu n = (2n+ l)7r/3 _1 which arises from the anti-periodic boundary conditions for fermionic fields in the path 
integral. The nucleon mass matrix is given here by 

M N = g(a + i'Ysnn) (A4) 



for non- vanishing <f> a — (cr, tt). Inserting (A3) into (A2) and performing the determinant over Lorentz and isospin 
indices we find 

^baryon^) = -2 • 2 J d 4 kln[k 2 + g 2 ^] (A5) 

where fc 2 = u> 2 + k 2 . In this result, one factor 2 has arisen from determinant in the Lorentz indices and may be 
attributed to the spin degrees of freedom, while another factor 2 is due to the isospin degrees of freedom. Note that 
fermionic loop contribution is larger than that of bosonic loop by another factor 2 due to the distinction of particle 
and anti-particle in case of fermion. In the presence of non- vanishing pion field g 2 4> 2 plays the role of the square of 
the effective nucleon mass, as expected simply from the symmetry consideration. 

Adding this term to the effective potential gives rise to new terms in the stationary conditions for the meson fields 
and the meson mass equations. We could modify our "Hartree approximation" by inserting the following term to the 
meson self-energy to the right hand side of (|57j) and ([58]) : 

H^(m N )=Il^(m N ) = -8gH (m 2 N ) (A6) 



with 

1 



k 2 + M 2 



where the factor 8 = 2 • 2 • 2 in (A6) accounts for the particle-antiparticle, spin, and isospin degeneracy and the 
subscript F in (A7) indicates that the sum is taken over the fermionic Matsubara frequency. 

This result resembles remarkablly the ones we have obtained for the meson loop contribution to the meson self- 
energies (p9|). Inclusion of the baryon loop contribution to the meson self-energy is not a consistent procedure, 
however, in view of the procedure we have adopted to drop the diagrams in Fig. ||, since both diagrams generate 
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non-trivial momentum dependence or dispersion in the self-energy of mesons and thus again invalidates our ansatz 
( |56| ) for the meson propagators. 




B 



o 



(a) 



(b) 



FIG. 9. Baryon loop diagrams which contribute to the self-energy of sigma (a) and pion (b) 



APPENDIX B: HIGH TEMPERATURE EXPANSION OF ONE-LOOP INTEGRALS 



Phase space integral in one-loop diagrams contain definite integrals of the following form: 



l dx 



1 



^x 2 + n 2 e V x2 +^ 2 ± 1 
For example, the boson (fermion) loop contributions to the meson self-energies is 



d 3 k 1 



1 



1 



(2tt) 3 E k eP E * ± 1 r 2ir 2 ± 
while the pressure of ideal gas of boson (fermion) with mass M is 

d 3 k 



P(/3;M) = ±/3- 1 



(2tt) 3 



In 



lie 



-/3Vfc 2 +M 2 



(4) 



(MP), 



where +(— ) sign is for fermion (boson). 

The functions (p) satisfy the following recursion relations: 



-I 



(n-2) 



and have the following limits for fi — > 0: 



li n) (0) = r(n)C(n), 
/i n) (0) = (1-2 1 -™ 



)r(n)C(n), 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 
(B6) 



where T(x) is the gamma function and C( x ) = Yl^—i l/ mX i s the Riemann £ function. 

We wish to evaluate the integral I±\[i) for the small nonvanishing value of /j = M/3, corresponding to high 

temperatures (T >> ill), however, the power series expansion of I± (p) would break down at fi n since the gamma 
function T(x) is singular at negative integer values of x. These singularities originate from the infrared (small x) part 

(n) 

of the integral I± . In the case of boson it is enhanced due to the additional singular behavior of the distribution 
function of massless boson at small x: l/(e x — 1) ~ 1/x. 

for even integer value of n. Here we quote 

-1/2 



Series expansion of (p) was derived by Dolan and Jackiw in 
some of their useful results: 



r(0) 



7T 
%1 



1 n 1 1^1 



47T 



7T 1 U 1 

= 2^ + 2 ln 4^ + 2 7 



2 ^ n 

n=l 

03) , 

167T ^ 



47r 2 n 2 



1 



305) 4 
64tt 2 ^ 



o(m 6 ), 



(B7) 
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r(0) 



2 7T 2 



1 ^ 1 
7+0E- 



2 " n 

n—1 





-1/2 


+ 47r 2 n 2 ) 


- 1 





1 M 1 

2 7T 2 ' 



C(3) 2 3C(5) 4 



16tt' 



64tt 2 



/i 4 + 0(//>) 



(B8) 



where 7 = 0.57721 • • • is Euler's number and the numerical values of the £ function at relevant points are £(2) = 7r 2 /6, 
C(3) = 1.2020- • •, C(4) = tt 4 /90, C(5) = 1.0369- • ■, and so on. The series expansion of I { ± ] {p) and I^jji) can be 
obtained from (B7) and (BS) by integrating the differential equation (B4) with the boundary conditions (B5) or QB6|): 



The result is 



for boson and 



7fM=7^(0)-i^ 2 rf M ' 2 /£ 0) (M') 



' 2 1 2 V) 



7T~ 
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1 2, /» 
— u — -u In — 

2 P 4 P 4tt 



1 1 



7 A* 



2 C(3), , 



32tt 



M 4 + 0(M b ), 



3 4 M 
— u In — 
16 P 4tt 



7 U 4 + 0(M 6 ), 



r(2) 
! + 

r( 4 ) 



-2+7]^ 



2 C(3 V 4 + 0(M 6 ), 



32tt 



7^ 
120 



3 4 /i 3 

T^M in _ 

16 7r 32 



2 7 U 4 + 0( M b ), 



(B9) 
(B10) 

(Bll) 
(B12) 



for fermion. 



[1] S. Weinberg, in Chiral Dynamics: Theory and Experiments eds. A. Bernstein and B. Holstein (Springer, 1995) p. 3; H. 
Leutwyler, ibid, p. 14. 

[2] Y. Nambu, Phys. Rev. Lett. 4, 380 (1960); Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961). 
[3] J. Goldstone, Nuovo Cimento 19, 154 (1961); J. Goldstone, A. Salam, and S. Weinberg, Phys. Rev. 127, 965 (1962). 
[4] For example, see H. E. Stanley, Introduction to Phase Transitions and Critical Phenomena, Oxford University Press, Inc. 
(New York, 1971). 

[5] D. A. Kirzhnits and A. D. Linde, Phys. Lett. 42 B, 471 (1972); Ann. Phys. 101, 195 (1976). 

[6] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 

[7] S. Weinberg, Phys. Rev. D 9, 3357 (1974). 

[8] T. D. Lee and G. C. Wick, Phys. Rev. D 9, 2291 (1974). 

[9] G. Baym and G. Grinstein, Phys. Rev. D 15, 2897 (1977). 
[10] M. Gell-Mann, and M. Levy, Nuovo Cimento 16, 705 (1960). 
[11] B. Lee, Chiral Dynamics (Gordon and Breach, 1970). 

[12] For review of recent works, see: A. Ukawa, Lattice 89, Nucl. Phys. B (Proc. Suppl.) 17 (1990) 118; S. Gottlieb, Lattice 90, 
Nucl. Phys. B (Proc. Suppl.) 20 (1991) 118; F. Karsch, Lattice 93, Nucl. Phys. B (Proc. Suppl.) 34 (1994) 63. K. Kanaya, 
to appear in the proceedings of Lattice 95. 

[13] T. Hatsuda and T. Kunihiro, Phy. Rev. Lett. 55, 158 (1985); Prog. Theor. Phy. 74, 765 (1985). 

[14] S. Weinberg, Phys. Rev. Lett. 18, 188 (1967); Physics 96A, 327 (1979). 

[15] J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987); Nucl. Phys. B 307, 763 (1988); ibid 321, 387 (1989). 
[16] J.D. Bjorken, Int. J. Mod. Phys. A 7 (1987) 4189; J. D. Bjorken, K. Kowalski and C. Taylor, SLAC preprint, SLAC-PUB- 
6109 (1993) 

[17] K. Rajagopal and F. Wilczek, Nucl. Phys. B404 (1993) 577; S. Gavin, A. Gocksch, and R.D. Pisarski, Phys. Rev. Lett. 
72 2143 (1994); S. Gavin, and B. Miiller, Phys. Lett. B329, 486 (1994); M. Asakawa, Z. Huang, X. -N. Wang, Phys. Rev. 
Lett. 74, 3126 (1995); J. Randrup, Phys. Rev. Lett. 77, 1226 (1996). 



20 



[18] D. Boyanovsky, H. J. de Vega, and R. Holman, Phys. Rev. D 51, 734 (1995). 
[19] F. Cooper, Y. Kluger, E. Mottola, and J. P. Paz; Phys. Rev. D 51, 2377 (1995). 

[20] We follow the notation in C. Itzykson and J.-B. Zuber, Quantum Field Theory (McGraw-Hill, New York, 1980). 
[21] B. D. Serot, and J. D. Walecka, Advances in Nuclear Physics 16, eds. J. W. Negele and E. Vogt, Plenum Press, New York 
(1986). 

[22] A. Fetter, and J. D. Walecka, Quantum Theory of Many-Particle Systems (McGraw-Hill, New York, 1971). 
[23] J. Kapusta, Finite- Temperature Field Theory (Cambridge University Press, Cambridge, 1989). 
[24] R. Jackiw, Phys. Rev. D 9, 1686 (1974). 

[25] S. Coleman, R. Jackiw, and H. D. Politzer, Phys. Rev. D 10, 2491 (1974); M. Kobayashi, and T. Kugo, Prog. Theor. 

Phys. 54, 1537 (1975); see also, A. Kerman and D. Vautherin, Ann. Phys. 192, 408 (1989). 
[26] Ref. U appendix C. 

[27] J. M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. D 10, 2428 (1974); G. Amelino-Camelia and S.-Y. Pi, Phys. 

Rev. D 47, 2356 (1993). 
[28] H.-S. Roh, J. Korean Phys. Soc. (Proc. Suppl.) 29, S372 (1996). 



A m a (MeV) e = e/(mg<fo) 

50 404 0.219 

100 555 0.0846 

200 772 0.0378 

500 1208 0.0142 

1000 1703 0.0069 



TABLE I. A dependence of the sigma meson mass at zero 
temperature and the scaled symmetry breaking parameter 
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